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I Advantages of using a low-energy effective theory to study bound state properties are briefly 

discussed, and a nonperturbative implementation of such an effective theory is described within 
the context of nonrelativistic quantum mechanics. The hydrogen atom, in the approximation of 
a structureless, infinite-mass nucleus, but with the leading relativistic and radiative corrections 
' included, is used to demonstrate the construction and solution of the effective theory. The re- 

, suiting Hamiltonian incorporates a finite ultraviolet cutoff and can be solved nonperturbatively. 

' An appendix lists explicit formulae for the various matrix elements necessary to diagonalize the 

■ Hamiltonian using gaussian basis sets. 

o 
o 
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1 Introduction 

^ , The study of QED bound states tests our understanding of nonperturbative field theory, and of the 
QJ , techniques used to study these nonperturbative phenomena. Comparison of theoretical predictions 
to precision experiments gives the best determination of various fundamental constants (e.g. Roo 
from hydrogen spectral lines, from the muonium hyperfine splitting). Once all of the relevant 
^ ■ constants are determined, a disagreement in the theory-experiment comparison can act as a signal to 
■ new physics. 

Given these motivations to study precision bound state QED, there are some difficulties which 
must be addressed. Two of these difficulties are the complexity of relativistic formalisms, and the 
cumbersome nature of bound state perturbation theory. 

Traditionally (although this tradition appears to be shifting), precision bound state calcula- 
tions have been performed in relativistic formalisms, using the Bethe-Salpeter equation or one of its 
variants. A drawback of the integrated approach of such a formalism, which does not distinguish 
between low-momentum and high-momentum contributions (since this would violate relativistic in- 
variance), is the difficulty in choosing an optimal gauge. For instance, Coulomb gauge could be most 
effective for parts of the calculation involving low momentum, but is cumbersome for high-momentum 
where Feynman gauge may be more appropriate. However in Feynman gauge, spurious low-order 
contributions appear which, though cancelling in the end, make it difficult to determine exactly which 
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Feynman diagrams must be evaluated to obtain a given accuracy. When a low energy effective theory 
is used to separate high- and low-momentum modes in a systematic way, the appropriate gauge can be 
chosen independently for both momentum regions. Technical simplifications of this sort result from 
the effective theory taking advantage of the essentially nonrelativistic character of the system. 

Bound state perturbation theory involves nontrivial sums over intermediate states, and in high 
orders the terms to be evaluated increase in number and complexity. Another problem is the reliance 
on an unperturbed state which forms the basis for perturbations. This is not such a problem in the 
single particle case, where the familiar Schrodinger-Coulomb wavefunctions are an obvious choice, 
but for multiparticle systems such as helium no analytic unperturbed wavefunction is known. These 
difficulties can be overcome by solving the effective theory nonperturbatively Q 

In the rest of the talk, after mentioning some examples of low-energy effective theories, I 
demonstrate, in stages, the construction of an effective Hamiltonian for the hydrogen atom with 
leading relativistic and radiative corrections, in the approximation of a structureless, infinite mass 
nucleus^. Ref. gives a more modern application, to the decay rate of orthopositronium. 



2 Low-Energy Effective Theories 

An example of a low-energy effective theory is NRQED (nonrelativistic QED) field theory, described 
by the Lagrangian 

Aff = -\{E^-B^)+i:^(^dt-e4> + ^ (1) 

+ (^ ^ x) 

d d 

together with a cutoff prescription, at momentum scale A. Here = V -|- ieA is the covariant 
derivative and V; X Pauli spinor fields, ci, C2, . . . , di, (i2i • • • are renormalization constants which 
must be determined by matching the predictions of NRQED to those of QED. The local operators 
parameterized by di, ^2, . . . account for the short-distance/ high- momentum states which are excluded 
by the cutoff |]. 

Another, and familiar, example of a low-energy effective theory is nonrelativistic quantum 
mechanics, described by a Hamiltonian, 



r, ~r Vlong— ^rangc 

2m 



HcS = £- + ^long^ra„gc(r) + -4,5i(r) + -^(-V^5i(r)) + -^p.5i(r)p + .... (2) 



^long-range includcs all long range interactions (the analogue of those operators parameterized by 
Cj in Eq.(l^) ). The remaining operators are the leading terms in an expansion of local operators. The 
structure of V|ong-range(^') and the values of the coefficients di are determined by matching predictions 

^ The low-energy effective theory will still be determined perturbatively, as an expansion in a and the typical atomic 
velocity v. It is the effective theory, and not QED itself, which will then be solved nonperturbatively. 

^ The ideas here are in most cases a rewording, and in others a slight extension of the presentation in Ref.||l|. 
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of the effective theory to those of the true theory (in the present case, QED) 0. It is on the effective 
Hamiltonian which we will focus, and for which we will formulate our nonperturbative implementation. 



3 The Schrodinger-Coulomb Problem 

First, consider the Schrodinger-Coulomb problem, where the "true" theory has Hamiltonian 

^sc^f-<^. (3) 

Zm r 

We will treat this as an example problem, to demonstrate the basic ideas in a familiar setting. 

Our goal here is to build an effective Hamiltonian in the form of Eq.(|2|), which describes 
the theory defined by Hsc- We first introduce a cutoff A ~ m which effectively removes states of 
momentum p > A from the theory. This is reasonable, since for such large momenta, e'^e" pair 
creation and other relativistic effects become important in QED 0. The particular cutoff we choose is 
in the form of a gaussian multiplying the momentum space potential: 



1 /I 



Ait 47r 

where the second line is the Fourier transform of the first: 



(4) 



1\ r d^a • 47r 9^ 1 Ar 
\ / " y e^q-r ' ^ -erf(^), (5) 



with erf (x) the error function. Introducing the cutoff modifies the Coulomb potential at small distances 
r <^ 1/A. With this form of the cutoff, a natural choice for the local operator 5\{r) is: 

Si (r) - [ ^e^q-e-i - (q) 

The other local operators in Eq.(|^ are derived from 5\{r) by differentiation. 

With our modified Hamiltonian, we no longer expect to solve the eigenvalue problem for the 
energies analytically; instead, we look for a numerical solution Q Two approaches can be implemented 
easily: 

(i) Differential Equation Integration: With i^gff expressed as a differential operator, the 
Schrodinger equation is: 

(^-^V^ + V{r)^i;{r)=Ei;{r). (7) 

Together with appropriate boundary conditions {ip{r — > oo) = 0, ip{r — > 0) finite), this eigenvalue 
problem can be solved using standard numerical differential equation integration routines. 

^ Alternatively, NRQED field theory could be used in an intermediate step: NRQED matched to QED, QM matched 
to NRQED. 

* It happens that Hsc is well defined without a UV cutofi^, so that introducing one may seem artificial here. However, 
when relativistic corrections are included, a cutoff is essential. And even in the present case, there are numerical 
advantages to having the cutoff in place — for example, fewer basis functions are required in the matrix diagonalization. 

^Since an analytic solution takes the form of a series expansion in a, a numerical approach is more natural as a 
nonperturbative tool. 
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(ii) Matrix Diagonalization: This technique will be our main focus, it being more readily adapt- 
able than the differential equation method when radiative corrections are introduced. Elements of a 
(linearly independent, but not necessarily orthogonal) basis are first chosen, for example: 



^^{v) = YUe,ct>ye ^-r. (8) 

where is a spherical harmonic and Ri runs over an appropriate range of distance scales. We now 
take matrix elements of all operators involved: (for clarity, we replace the three indices i,l,m above 
by a single index) 

Now taking matrix elements of the Schrodinger equation, we obtain: 

^ {p'^y^ + v'A = Ew'^ijj , (10) 



2m 

or as matrices: 

( 

. 2m 



— p2 + y ) V = EWi;, (11) 



where ^ = ^■'ipj, and repeated indices are summed over. Since we have not taken our basis elements 
to be orthonormal, W^^ ^ S^^ . To recover a more standard form, define, as matrices, 

W'Ww^^^ = V 

W2tp = tp. 

Then Eq. (|Tl|) , upon multiplying both sides by W^2 becomes: 

1 



p'^ + V]Tp = Eip. (13) 
2m / 

This matrix eigenvalue equation can be solved easily using standard matrix diagonalization routines^. 
The long range potential of Eq. (^) is 

^long-range (0 = -{Zo) (^^)^ = Vc- (14) 

To complete the determination of Heg, we must evaluate the necessary local operator coefficients 
di, (^2, . . .. First, we should decide which coefficients must be evaluated, and to what accuracy. This 
is accomplished by noting that 

{^o\^Si{r)\^Po) « 4l^oP « d,ma\ (15) 
m^ m^ 



See for example Ref. Q . For the relatively small computer time required for the problems treated here, a commercial 
program like "Maple" or "Mathematica" works fine. 
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So, to determine energy levels through, say, O(ma^), we need determme di through O(a^). Terms 
parameterized by ^2 and d^ contain two more powers of p/m « a, and so contribute at 0{d2ma^), 
0{d3ma^) to the energy. 

We will implement the matching procedure between the true and effective theory by perturba- 
tive matching of scattering amplitudes for the process k ^ 1. The lowest order scattering amplitude 
(the Born term) is simply the momentum space potential: (here q = 1 — k is the momentum transfer) 

d^ o d^ 



(i|r|k) 



AniZa) di ao o 

+ ^ + + — 



47r(Za 



Ik 



<3 



(16) 



+ 



di 2Tr{Za)\ f d2 



+ 



H 



7r{Za) \ 
2A4 ) 



The first term matches the 0{a) result of the "true" Coulomb theory. Requiring the other terms to 
vanish, through 0{a), gives: 



7(1) 



-27r 



m 
A2 

4 

vr m 
2 A4 



(Za)4^^ + {Zafd^2^ + (Zafd^^^ + 



{Za)d^2^ + . , 



where di _ 

At 0(0^), our power counting shows that we need only evaluate d^^ 
threshold scattering amplitudes in the true and effective theories |^ 



(17) 

(18) 
(19) 

(Za)4^^ + •••• 

. To do so, we consider 




Requiring this quantity to vanish yields: 



7(2) 



vr 



10 

y 




47r m 



96 V A 



By a similar calculation at 0{a^), with the above calculated values of d^^ and d[^^ 0. 



7(2) 



d 



(3) 



'bVS 1633 




(20) 



(21) 



(22) 



A comparison of the first few S -state energy levels calculated both with and without counterterms 
are shown in Table |l|. As expected, with counterterms included the energy levels are in agreement with 
the Coulomb spectrum through 0{ma^). This matching procedure can be extended systematically to 
higher orders. 

^ Here we are using the relation Tik — Vik + VipGp{E)Tpk, where Gp{E) = {E ~ /{2m) + ie)^^ is the nonrelativistic 
propogator and for the scattering state, E = fc^/(2m). 

* For energy levels correct through 0{ma^), it would be sufficient to include only cf^^ and d^^\ with d2 — 0. The 
value of d^i'' would then be given by Eg. (pi)) at = 0. 
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n 


E/Eq 


E/Eo 


1 


0.99850 


1.0000000035 


2 


0.24981 


0.2499999996 


3 


0.11106 


0.1111111109 



Table 1: S -state energy levels at a = 0.02, A = m. Here Eq = —mo? 12 is the ground state Coulomb 
energy. In the first column, d\ = di = 0. In the second, these parameters take the values calculated 
in the text. 



4 Relativistic Corrections 

The leading relativistic corrections to the Schrodinger-Coulomb Hamiltonian take the form [||: 

I ^^^"^ j3(r) I (^")^"'^ (23) 
'^rrfi 2m^ Am? r'^ 

We can immediately see problems with these corrections as they stand, if we are to solve 

nonperturbatively. At high momentum, the p"^ term dominates over p^, causing the Hamiltonian 

spectrum to be unbounded from below. Also, the (5-function is too singular and all second-and higher 

order perturbations involving this term will be divergent. To remedy these problems, we introduce a 

cutoff: 

S^r) ^ 6l{r), (24) 

1 1 /i\' 1 ^l\ , , 

(25) 



giving the cutoff Darwin and spin-orbit potentials: 

4m^ r \r J \ 

For the term, we choose to work with an equivalent form |^: 

/ ^ {2mf{E - Vf, (28) 
which defines our kinetic energy correction potential: 

Vk{E) = -^{E-Vc)\ (29) 

where Vc is the cutoff Coulomb potential. 

Working through 0{ma^), the only counterterms necessary are which is given in Eq.(fi!7D, 

(2) 

and d\ , which must be recalculated with the relativistic corrections in place: 

12) ^( 10 {m\^ m 1 A\ 



^Here we could also apply the cutoff directly to the operator, for example p*exp(— p^/A^). The other form happens 
to make evaluation of the counterterms more convenient since then only the free propogator /{2m) appears in scattering 
calculations. 
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(This calculation is similar to that of Eq.(|20|), but with T^rue now referring to the scattering amplitude 
of QED in the external field approximation.) 



5 Radiative Corrections 

It remains to include the effects of radiative corrections — in the language of Coulomb gauge, the effect 
of the transverse radiation field. Our analysis is complicated here by the fact that soft photons, of 
energy E ~ ma^, can no longer be described by instantaneous potentials in the electron's Hamiltonian, 
since their characteristic propogation time {6t ^ 1/E) can be comparable to bound state timescales 
{6t ~ l/(ma^)). We can deal with this by expanding the state space to include two channels: one 
with just the electron, another with the electron and a transverse photon. The wavefunction now has 
components in both sectors: 

l-^it']- (31) 

The new terms in the Hamiltonian which describe the coupling between channels are: 

6H = H^-—p- A(r), (32) 
m 

where 

^7 = / 7^ E '^'i «!l,e(q) «q,^(q) (33) 



e(q) 



is the photon kinetic energy operator, and 

AW = /|fji;;^(W- + '-)^ (34) 

For photons of momentum q <C ma (i.e. q < ma^), the multipole expansion is valid: 

A(r) = A(0)+r- VA(0) + ... (35) 

For q » ma^ (i.e. q ^ ma), the photon propogation is effectively instantaneous, and can be described 
by instantaneous potentials in the electron Hamiltonian. So, denoting the potentials already present 
(Coulomb and relativistic corrections) by V, our coupled channel Hamiltonian takes the form: 

2 

H^^ + V + H^-—p- A{0) + ... (36) 
2m m 

where the dots represent any new instantaneous potentials generated by interaction with the transverse 
photon sector. 

The coupled channel problem can be reduced to an effective energy-dependent Hamiltonian 
acting only on the electron subspace. In general, if the coupled Schrodinger equation has the form: 

( f't H^+H, ) ( ( ) ' ^^^^ 

then the effective Hamiltonian is 

H = He + H'{E - He- H^y^H'^. (38) 
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In the present case, 



He 

H' 



p 

2m 

e 

ni 



+ V 



p- A(0), 



(39) 



and is given by Eg. (|33D . Eq.(^) can be evaluated explicity. Including radiative photon modes of 
momentum ^ g < ~ m, and assigning the photon a mass A to control infrared divergences, 



H 



a 



(40) 



+ 



^1 



E-H, 



+ 1 



with qx = \/ q'^ + A^. 

The final step in our analysis is to compute any remaining instantaneous potentials, by compar- 
ing our effective theory to QED. We do this by examining the 0{a) correction to the QED scattering 
amplitude, which can be written in terms of vertex form factors, and a vacuum polarization contribu- 
tion: 



t true 



-AniZa) 47r(Za) 



with 



F2 



Am? 



a 



1 



2i 1 X k • fj 



Fa + 



-47r(Za) 



a 



IStt m^ 



a 
2^ 



In 



m 



+ 0{\)]+0{q' 



To this we must compare the 0{a) contribution coming from Eq.(|40|) P^ : 

-47r(Za) 



Q 9 

q 



Comparing Eqs.(|4l|) and ([44D, we see that potentials: 

4:a{Za 



Vfi 

Vf2 

VvP 



{Za) a 
4m2 27r 
4a (Za 



m 11 \ , , 

In — — H \5l (r) 

2A^ 24 / ^ 



4vr<5i(r) 



2I riiL.. 



(41) 

(42) 
(43) 

(44) 

(45) 

(46) 
(47) 



associated with Fi and F2 form factors, and vacuum polarization, respectively, must be added to our 
effective theory. Having determined these potentials, we now set A = 0, and evaluate the 0{a) piece 



instead of « ma avoids the appearance of factors In a in higher order contact terms, but is not 



^° Taking 

essential. The difference will be accounted for by the remaining instantaneous interactions — see Eq.(|45|). 

Here wavefunction renormalization must be taken into account, since the potential is energy-dependent. 
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a 


(S(l, 1/2,0) -£;'(!))/(- 


-4a5/37rn3) 


(-E(2, 1/2,0) -^'(2))/(-4aV37rn3) 


0.04 


2.79656 




2.51928 


0.02 


2.88486 




2.65806 


0.01 


2.92986 




2.72998 



Table 2: S -state energy levels with relativistic and radiative corrections, at A = = m. Here 
E'{n)lm = -{ZG?)l2n^ - {Zaf/2n^{l - 3/4n) + a(Za)^/7rn3(8/3 In a'^ + 10/9 - 4/15) is the part 
of Eq.([50|) excluding In/co- At a ^ the second and third columns converge to lnA;o(l,0) = 2.98413 
and ln/co(2,0) = 2.81177, respectively. 



ofEq.® B 

Vls{E) = ^e-^ f {^ + yc- E) In 1 p'e-^ . (48) 

37rm"' 2m p_ _|_ T/r^ _ p 

We can now present the complete effective Hamiltonian, including leading relativistic and 
radiative corrections: 

2 

HeE = -^ + Vc + Kel + ^rad + , (49) 

where Vc is the cutoff Coulomb potential, l^ei = ^A'(-E') + Vd + Vso incorporates relativistic 
corrections, T4ad = ^LS'(^) + Vpi + + Vvp contains the radiative corrections, and Vet = 
{(f^^ /m? + {Za)d^^^ /rn?')5\{r) is the counterterm potential. For this application, we have (Eqs. 

(0),® ) 



n 

.(2) 



A2 



/ 10 (my m 1 A 



The spectrum of this Hamiltonian can be readily determined using matrix diagonalization |^, 
and reproduces (Table ^) the well-known result through O(ma^) Q : 

^ ■' ' 2n2 2n3 Vj + 1/2 ^n) ^ ' 

ma{Zaf /8, 10 4 

+ 5^,0 -In Za i + — - — 

'nrr V \o 9 15 

+ (1 - (^io) — — 777 , , lnA:o(n,0 

" '"^ V 2^(^ + l)(2/ + l) J 3 ' ^. 

It is worth noting that higher order terms , including ones containing ln(a), appear automat- 
ically when the theory is solved nonperturbatively. This can be important in high orders when the 
appearance of factors (In a)", for large enough n, causes poor convergence of series expansions. 

The first term in square brackets in Eq. ( |4o| ) is absorbed by a mass renormalization. "LS" stands for "Lamb shift", 
since this potential accounts for the dominant part of this effect. Note also that a cutoff in the form exp(— p^/2A'^) has 
been included in Vls- 

The necessary matrix elements of the various operators between gaussian basis functions are listed in the appendix. 
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Appendix: Matrix Elements for Gaussian Basis 

Following is a list of matrix elements between basis functions 



(51) 



It is convenient to define the quantity Rij = (R^ ^ + Rj ^) Then the matrix elements are: (a 

factor Smi,m.jSi^,L is suppressed) 



p 

i) 

Slir) 



^ + 3/2 r 

1 



1 / 1 



+ 



2\Rf i?f 



A2(«+i) ■ 



da 



a=l/(A2K2^.) 



2m/2 12-^.3 

-7TT\ -T { -+l 

47rA2' V TT V2 



1 + 



WW. 



1+3/2 



3A2 



5l(r) 



2'+3/2 /2 /5 , 
-r - + / 



47rA2^-2V7r \2 



1 + 



1+5/2 



(52) 
(53) 
(54) 

(55) 
(56) 



10 



(58) 



For the operator appearing in VlSj we concentrate on the ground state, with / = 0. To first 
order in Vls, we need only consider couphng between 5- and P- states {S ^ P ^ S transitions)]^. 
Here the sum over spherical harmonics results in the replacement J2iP^i''')P^ ~^ ~(^Pr)(" ■ OC^Pr-); 
where ipr = d/dr, and: 

= ^,1^,0 + . (59) 

Radially excited states (/ > 0) can be treated similarly; for instance / = 1 requires a sum over couplings 
between P- and S- and between P- and D-states. 

The Ri should take values which cover the range from ~ 1/A to ~ l/{ma). For example, 
the values of Table |l| were generated using Ri = Rmm{Rma,x/ Rmin 

yl'^ , with i = {)..N, N = 60, 

-Rmin = 0.1/A, i?max = 100/(ma) and A = m, a = 0.02. For problems with more than one angular 
momentum channel, a separate set of basis functions should be used for both channels, e.g. Ng S-states 
and Np P-states for the hydrogen atom problem with radiative corrections. 



Higher orders will introduce S^P^D—*P—^S transitions, but these are suppressed by several powers of a. 
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